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The problem considered is that of determining the failure statistics for a system 
consisting of a single part and n spares, in which it is assumed that that failure occurs at 
different rates in use and in storage. The resulting equations can be solved explicitly 
when both types of failures follow a Poisson law. Otherwise algorithms are given for the 
calculation of the failure probabilities. 



The problem which we consider in the following 
paper is the calculation of the lifetime statistics of a 
component with n spares with the proviso that re- 
dundant elements can fail in storage. It is assumed 
that failed components are instantaneously replaced. 
Similar models have been considered before in the 
literature, particularly by Proschan and collaborators 
[1,2] \ who were interested in determining the opti- 
mal spare parts kit for a system of components in 
series, given a fixed amount of capital to provide the 
spare parts. If it is assumed that components do not 
deteriorate in storage then the calculation of tlie dis- 
tribution of time to failure is simple since it is just 
the sum of individual failure times. In the present 
situation tlie calculations are somewhat more com- 
plicated, and we have only succeeded in finding a 
solution when the reliability functions for on-line 
failures and for in-storage failures are exponential. 
In other cases we prescribe algorithms for the 
solution. 

The first case to be treated is one in which it is 
assumed that there are two failure rates independent 
of time. Let fxidt be the probability that a com- 
ponent will fail in (t, t+dt) in the stockpile and fx^dt 
be the probability that a component in service will 
fail. Further, let 

Pj{t)=Pr{j components in the system at 

time t of which one is in service}. 

This is a Markov process, for which the state equa- 
tions are easilv seen to be 



1 Figures in brackets indicate the literature references at the end of this paper. 



PoW=M2Pi(0 



P,{t) = (/:m. +M2)P*+, (0 -{(k- 1)m. +M2)P.«) 



P.W = -((«-i)mi+M2)P„(0 



(1) 



in winch it is assumed that there are initially ri com- 
ponents in the system, one in use and n— 1 in storage. 
The state equations are most easily solved by intro- 
ducing Laplace transforms: 



^.(5)-J^ e-''. 



PS)dt, 



(2) 



For simplicity we will also introduce the notation 

o-m=^Mi+M2. (3) 

A recurrence relation is easily written down for 
the ifkis). It is 
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<Pm + l{s) l<m<7l 



(4) 
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Let us denote the probability that the system is 
down to k components for the first time in {t, t-\-dt) 
by Trfc{t)dt. Then it is clear that 



'jrS) = a,P,+,(t) 



or if 



then 



<(s)- 
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''Mt)dt 



7r?(s)=V 



J=k S + Cj 



(5) 



(6) 



(7) 



In particular the Laplace transform of the probability 
density for the time to system failure is just ^^{s). 
The Laplace transform of eq (1) is easily inverted by 
expanding it into partial fractions and inverting the 
resulting series term b}^ term. We find that 



(n—l \ n — 1 /n—1 
n <r,)2](n' 
j=k / r = k \m=k 



1 



1 



(o-^— 0-,)/ (s + ar) 



(8) 



where the prime denotes the fact that the term 7n=r 
is omitted. The inversion of 7r*(5) yields the follow- 
ing expression for x^(0 



TTjcit) =a^2 



where 



(,,l)(.^,,l).,. (,_.+!) 



(n-l-k)\ 




g -^2 (l+ta)«(l — e -ap-it^ n- 


■-* (9) 


a=Mi/M2. 


(10) 



In particular the density for the time to SA^stem 
failure is 



To{t)-- 



'M2 / ^ -.N. 



(7^-1)! 



e-t.,t(l-e-af.,t)n-l (11) 



which can also be derived more directly. Notice that 
in the limit a-^0 this expression goes over into 



Toity 



>M2 



(n-l)\ 



e~f^2t 



(12) 



which is valid when there is no deterioration in 
stockpile. Nevertheless it is valuable to have the 
expression of eq (8) for the Laplace transform of 
TTjclt) since it can be used as a moment generating 
function. 

Moments of the time to failure are found through 
the relation 



r==(- 



-ir^^^mu^ 



o+- 



(13) 



For the mean and variance of the time to failure we 
find 



7 ' 



n-l 



1 



^i2 7=0 1+^Q: 
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1 n—l 



AU (H-i«)' 



(14) 



When a is equal to zero, which is equivalent to saying 
that there is no deterioration in storage, t and a^ are 
proportional to the original number of components. 
When a is finite and n increases, the approximate 
values of t and o-^ are 



?^-— log 71 -f 0(1) 
Ml 



Ml 



(15) 



where l{x, a) is the generalized zeta function [3]. 
The result of significance is that for finite a, t increases 
at a much slower rate with n than it does when there 
is no deterioration in storage. At the same time the 
variance remains bounded as the number of compo- 
nents increases. 

A plot of t(a)lt(0) is shown in figure 1 as a function 
of a for some representative values of n. This rep- 
resents the ratio of mean time to system failure with 
and without failure in storage. It is obvious that 
the larger the initial number of components, the 
more significant is the effect of failure in storage. 

Another result of some interest for this system is 
the reliability function E{i) : 



Bit) 



=j; 



TTo (7)^/7 



= M2 



r(n)r 1+ 



a/ 



g-M2^(l— g-«M2^)^-lrf^ 



=1 



exp(-aiU2 



-d"') 



where 



Ip(n, m) = 



j: 



x''-\i—xy''-kix 



C ^n~l 



(16) 



(17) 



(l-x)'^-'dx 



is the incomplete beta function. When n is large 
and a is small (both arguments in the incomplete 
beta function are large) one can use the normal 
approximation to the incomplete beta function to 
find the formula 



R(t)^eric 



a \ a/ 



^ — q:jU2^ 



UH) 



g -aH2i f\ — g""«A'2^) 



(18) 
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E(a)/E(o) as a function of < 



= 2, .3, 5. 



erfc 



1 C" 

-\I2t J X 



-''"'-hit. 



(19) 



For small n one can use Pearson ^s tables [4], to eval- 
uate the reliability function. Several reliability 
functions are plotted in figure 2 for different values 
of n and a. When a = 0, Hit) is given by the expresson 
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(20) 



It can be seen that the effect of even a small amount 
of deterioration in storage is quite pronounced in 
terms of system reliability. 

We have not succeeded in fuuliiig a general solu- 
tion to eq (1) when /xi and /xo are allowed to be func- 
tions of time. In the special case 



/^i(0 



^2{t) 



=q:= constant 



(21) 



one can relate the solution to the one obtained for jUi 
and 1x2 a constant. For this purpose the new time 
variable 



^x 



T= H2(x)dx 



(22) 



is defined. In terms of this variable the equations 
for the probabilities read 

f=A(o 



■{ik-l)a+])I\ir),k=2,'S, 



dP„ 



= -{(n-l)a+\)P„{r) 



(23) 
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Figure 2. Plot of R(t) for 11 = 1, a=0 

11-3, a = 0.2 

and the solution of eq (1) given in eq (11) remains 
valid except tliat t is replaced by r. In this situation 
the expressions in eq (14) for the moments are no 
longer correct because of the change in the time 
variable. Wlicn jlii(/) and 1^2(1) do not satisfy eq 
(21), the difference-differential equations can be put 
in a form whicli is more convenient for a recursive 
calculation. With the reduced time being defined 
by eq (22), the eq (23) remains valid provided that 
a{T) is regarded as a function of r. If we now 
define a new set of dependent variables Ukir) by 
the prescription 

?7*(r)=P.(r)e'+«-"/;<'<^"'^ (24) 

then the Uic(t) satisfy the initial conditions 

UM=PM=hn (25) 

and the equations become 

U,{r) = (kair) + l)e-/o«^^^^^^,+, (r) 

k=0, 1, 2, . . .,n-l (26) 

from whicli the U's and finally the P's can be 
determined in succession. 

When the reliability functions are not exponential 
the problem becomes much more difficult and we 
have not been able to resolve it in closed form. 
However it is possible to derive a set of recurrence 
relations from which the system reliability can be 
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found. Let us assume that the reliability of a 
component in stockpile is Rs(t) and that Rq{t, x) is 
the probability that a component placed in service 
at time r will last for a time x or greater. Let 



(t{t, x) = - 



Mio 



dx 



(t, x) 



(27) 



be the probability density for the failure of a compo- 
nent in (x+T, x+T+dx) given the installation of the 
component on-line at time r and let ps{t) = —dBs(t)/ 
dt. Finally define a function o)]c{t) by 

cofc(^)(/^=Pr{ Failure of the on-line component in 

{tj t^dt)\k components in stockpile} . 

Then the cofc(^) satisfy the recurrence relation: 

[iiJ[-IST"* <-' 



^kK 



which is derived by considering two possibilities. 
Either the on-line failure is the first such which 
occurred, leading to the first term in the right hand 
of this equaticn, or the last on-line failure occurred 
at time r and there were k^j components in storage 
of which k survive the additional time t — r. The 
first passage time density tor time to density failure 
is just 

7ro(0 = ^o(0. (29) 

The formula of eq (28) is almost useless from the 

computational point of view, so that no detailed 

investigations of its consequences have been 
undertaken. 
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